Energetics of nitrogen incorporation reaction in SiC>2 
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We study using first-principles calculations the energetics, structural and electronic properties of 
nitrogen incorporation in SiC>2. We consider NO, NH, N2 and atomic N as the nitriding species 
interacting with a Si-Si bond of an otherwise perfect Si02 network in order to simulate the nitrogen 
incorporation near Si/Si02 interface regions. We find that all the species react with the Si-Si bond 
forming bridge structures with the Si atoms without dissociating, where NH and atomic N form 
the most stable structures. Concerning the electronic properties, our results show that incorporated 
NH is the only structure which does not introduce trapping center at the interface. The structures 
involving NO and atomic N are acceptors, whereas that involving N2 may be either a donor or an 
acceptor. The hydrogen passivation of the electrically active centers is also discussed. 

PACS numbers: 61.72.Ji, 71.55.Ht, 71.15.Mb 



Oxynitridc films have been extensively studied in the 
past years since they improve the reliability of mctal- 
oxide-semiconductors gate insulators. The main benefits 
of incorporating nitrogen into ultrathin SiC>2 films are 
the reduction of gate leakage currents and the resistance 
to boron penetration^]. The growth of ultrathin oxyni- 
tride films strongly depends on the reactant species {e.g., 
N 2 0, NO, NH 3 , N 2 ) and the technique used. Nitrogen 
can be incorporated into SiC>2 using either thermal ox- 
idation and annealing or chemical and physical deposi- 
tion methods. Thermal nitridation of SiC>2 in NO and 
N2O generally results in a relatively low N concentra- 
tion at the near-interface (Si/Si02) region^, The 
N incorporation is commonly associated to the reaction 
of NO molecule with Si-Si bonds at the interface, after 
diffusing through the oxidej^. On the other hand, N 
incorporation via annealing in NH3 is responsible for a 
relatively high N concentrations into the films 3] . This 
method provides both near-interface and near-surface ni- 
tridation, suggesting different nitriding species derived 
from NH3. More recently, higher N concentrations and 
controlled distributions have been attained by plasma as- 
sisted methods, typically using ions and radicals derived 
from N2 and NH3 as nitrogen sources 5] . Although the 
control of both the density and the distribution of nitro- 
gen into Si02 have been achieved at few-layer level, less 
is known about the early stages of the N-incorporation 
reactions at atomic level. 

In this work the energetics and structural properties 
of near-interface nitrogen incorporation are studied from 
first-principles total-energy calculations. We have con- 
sidered the N2, NH and NO molecules as well as atomic 
N as the precursor species, reacting with the Si-Si bonds 
in Si02 in order to simulate suboxide (Si 3+ ) or near- 
interface N incorporation. Our calculations were per- 
formed in the framework of the density functional theory, 
using a basis set of numerical atomic orbitals as imple- 
mented in the SIESTA codeQ. We have used a split- 
valence double-^ basis set plus the polarization functions 
as well as standard norm-conserving pseudopotentialsQ. 



For the exchange-correlation potential we adopt the gen- 
eralized gradient approximation^. We used a 72-atom 
ev-quartz supercell and the T point for the Brillouin zone 
sampling. The positions of all the atoms in the supercell 
were relaxed until all the force components were smaller 
than 0.05 eV/A. We also consider neutral and singly 
charged species, where the neutrality of the cell is always 
maintained by introducing a compensating background 
charge. Spin-polarization effects, which are important 
for the correct description of atomic and molecular re- 
action processes in SiO2 0], are included throughout the 
calculation. 

We study the chemical reactions occurring when the 
nitriding species NO, NH, N 2 and atomic N approach 
the Si-Si bond in the otherwise perfect Si02 network. 
The Si-Si bond in Si02 is formed when an O atom is 
removed from the network characterizing an oxygen va- 
cancy. This local geometry is close related to those found 
near the Si/Si02 interface^. Hereafter, the incorpo- 
rated species will be called [NO], [N 2 ], [NH] and [N]. For 
these reactions we only consider neutral species in their 
ground-state spin configurations. Initially, we study the 
structural properties of the nitriding species inside the 
largest interstitial site of a perfect Si02, exploring pos- 
sible reactions that they may undergo with the network. 
Our results show that none of the species considered in 
this work react with the network, remaining at the in- 
terstitial sites. This suggests that neutral species would 
be diffusing species in Si02- However, when we put the 
species close to the Si-Si bond, they are quickly incor- 
porated into the network forming the stable structures 
shown in Fig. ^ We now describe in details our results 
for each incorporated specie. 

The NO molecule shows two stable structures after re- 
acting with the Si-Si bond. In the lowest-energy struc- 
ture [NO], the N atom is threefold coordinated bonding 
with two Si atom, keeping the bond with the O atom, as 
shown in Fig.^a). Here, the Si-N and N-0 bond lengths 
are 1.78 and 1.32 A, respectively. The binding energy of 
NO, calculated as the difference in energy between the 
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FIG. 1: Local equilibrium geometries for neutral nitriding 
species after reacting with Si-Si bond in SiC>2. (a) the most 
stable geometry for the incorporated NO. (b) the second sta- 
ble geometry for the incorporated NO. (c), (d), and (e) are the 
stable geometries for the incorporated NH, N2 and atomic N, 
respectively, (f) the Si-Si bond or the suboxide region (Si 3+ ) 
in Si0 2 . 



interstitial and incorporated configurations, is found to 
be 3.8 eV. Whereas, the binding energy of the O atom 
in the [NO] structure is 4.2 eV. Thus, it would be eas- 
ier to remove the entire NO molecules than the single O 
atom. Fig.^b) shows the second stable geometry [NO]'. 
Here, NO is incorporated into the Si-Si bond forming a 
structure similar to the peroxyl bridge of oxygens PUpf. 
This structure is 0.65 eV higher in energy than the most 
stable one. The Si-N, N-O, and O-Si bond lengths are 
1.77, 1.40, and 1.71 A, respectively. The energy barrier 
for the NO molecule to change from the metastable struc- 
ture [NO]' to the most stable one [NO] is estimated to 
be 0.6 eV. The lowest-energy [NO] structure is an elec- 
trically active center, exhibiting a half-occupied energy 
level at 1.9 eV above the density- functional valence-band 
maximum (VBM) of a-quartz. The [NO] dangling bond 
may be passivated by capturing an H atom as recently 
suggested^. We find that the reaction involving the 
capture of an H atom from an interstitial H2 molecule 
is endothermic with an energy cost of 0.78 eV. However, 
the [NO] passivation may occur by the capture of an in- 
terstitial H atom. In this case, the reaction is highly 



exothermic with an energy gain of 4.1 eV. 

Figure^c) shows the equilibrium geometry for the NH 
after reacting with the Si-Si bond in Si02- We see that 
the N atom binds to both Si atoms forming the Si-N-Si 
structure, keeping the bond with the H atom. The Si-N 
and N-H bond lengths are 1.72 and 1.03 A, respectively. 
The binding energy of the NH molecule is found to be 
6.8 eV which indicates that this radical would be one of 
the most stable species for interface nitridation. We also 
found that the H atom is strongly bound to the N atom 
in [NH], having a binding energy of 4.6 eV. 

The equilibrium geometry for the incorporated N2 is 
shown in Fig. ^d) . The N2 binding energy is found to be 
1.2 eV, relatively small as compared with the other nitrid- 
ing species. This suggests that N 2 may be easily removed 
from the interface escaping from the Si02 film. Addition- 
ally, the [N2] structure introduces an empty midgap level 
at 2.7 eV above the VBM of a-quartz. 

For the N atom we also find a very stable structure af- 
ter reacting with the Si-Si bond [Fig. df)]. The binding 
energy of the N atom is found to be 6.5 eV. [N] shows a 
half-occupied energy level at 0.2 eV above the VBM of 
a-quartz. This dangling bond may be passivated by cap- 
turing an H atom from an interstitial H2 molecule which 
dissociates after the contact with the [N] structure. We 
find that this reaction is exothermic with an energy gain 
of 0.82 eV, where the resulting equilibrium structure is 
similar to the [NH] structure [see Fig. [Jc)] P ms an inter- 
stitial H atom. If the [N] passivation occurs by capturing 
an interstitial H atom, the reaction is highly exothermic 
with an energy gain of 5.0 eV. 

The relative stability of incorporated nitriding species 
in Si02 in thermodynamic equilibrium is obtained by 
comparing their formation energies. As we are simulat- 
ing the nitrogen incorporation at suboxide regions, i.e., 
close to the Si/Si02 interface, we consider the three most 
relevant charge states in order to address possible charge 
transfer between the incorporated species and the doped 
Si substrate. Therefore, the formation energies (Ef) are 
calculated as a function of the chemical potential of the 
nitriding species (fJ-x) and the electron chemical potential 
(Me) as 

^/([x] 9 ,/i e ) = E t {[x] q )-E t {Si0 2 )+Lio-m+q ve, (i) 

where E t ([X] q ) is the total energy of the incorporated 
specie X in the charge state q, Et(Si02) is the total en- 
ergy of perfect a-quartz and no the chemical potential 
of oxygen, /zx and are fixed at the values of their 
stable gas phases. Our results for the formation energies 
of [NO], [NH], [N 2 ], and [N] are shown in Fig.H In the 
figure, the symbols indicate the transition states or the 
thermodynamics levels, whereas the slopes of the lines 
(positive, zero and negative) indicate the corresponding 
charge states. Following previous calculat ions \ 1 ll 1 1 2j . the 
transition states were aligned with respect to the experi- 
mental band-gap edges of a-quartz, by matching the cal- 
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FIG. 2: Formation energies of the nitriding species after re- 
acting with the Si-Si bond in Si02, as a function of the elec- 
tron chemical potential (fi e ) in the energy gap of Si02- £v'° 2 
[ef lC>2 ] is the experimental valence-band maximum (VBM) 
[conduction-band minimum (CBM)] of a-quartz. The dashed 
lines indicate the density-functional VBM and CBM and the 
shaded area represents the experimental band gap of bulk Si. 
The Si midgap is chosen as the zero energy of /j, e . 



species for interface nitridation having a binding energy 
as high as 6.8 eV, where the N atom is typically threefold 
coordinated. NO would also be a very stable species at 
the interface having a binding energy of 3.8 eV. Our re- 
sults agree well with the finding that NO diffuses through 
the oxide and reacts at the near-interface region where 
N is bound to two Si atoms and one O atomQ. How- 
ever, NO does not dissociate spontaneously after reacting 
with a Si-Si bond. N2 may also be incorporated at the 
interface, forming a stable structure with the Si-Si bond. 
However, it has a relatively low binding energy of 1.2 eV. 
This suggests that N2 may be easily removed from the in- 
terface. Atomic N does not react spontaneously with the 
Si02 network as commonly believed, suggesting that it 
would be a diffusing species in Si02 . The binding energy 
of the incorporated N is 6.5 eV, resulting in a very stable 
species for interface nitridation. Concerning the electrical 
properties, our results show that NH is the only species 
that does not introduce a trapping center after being in- 
corporated at the Si-Si bond. [NO] and [N] are acceptors 
whereas [N2] may be a donor or an acceptor, depending 
on the Fermi level position. However, [N] and [NO] may 
be passivated by subsequent hydrogenation processes. 
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culated (+/— ) transition state of the interstitial H atom 
in a-quartz with its measured value of 0.2 eV above the Si 
midg and by using the measured valence-band 

offset of the Si/Si0 2 interface of 4.3 eV0. In Fig. we 
also depict the experimental band gap of bulk Si where 
the Si midgap position is chosen to be the zero energy of 

Me- 

As we are computing the energetics of the incorporated 
species close to the Si/Si02 interface, the variation of jit e 
must be considered only in the energy range of bulk Si. 
Fig. [21 shows that negatively charged [NO] is the most 
stable structure at the interface. The (0/—) transition 
state is found at 4.3 eV above the experimental Si02 
valence-band maximum (e^ 2 ) . This suggests that [NO] 
might capture an electron from the Si substrate charac- 
terizing an acceptor center at the Si/Si02 interface. The 
same behavior is found for [N] which might be another 
acceptor center at the interface. For the case of [N2] we 
find both donor and acceptor character, where the (+/0) 
and (0/—) transition states are localized at £^ lC>2 +4.5 
and £^ lC>2 + 4.8 eV, respectively. The formation energy of 
[N2] is approximately 2 eV higher in energy than [NO], 
suggesting that [N 2 ] would be the less stable structure 
at the interface. Finally, [NH] is found to be the only 
electrically inactive structure at the interface, having a 
relatively low formation energy. This suggests that the 
NH molecule would be the best choice among the consid- 
ered nitriding species for interface nitridation. 

In summary, we find that NH would be the most stable 
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